arXiv:1504.02469vl [astro-ph.SR] 9 Apr 2015 


Astronomy & Astrophysics manuscript no. 25757 ©ESO 2016 

November 15, 2016 


Astrometric planet search around southern ultracool dwarfs 

Ml. Discovery of a brown dwarf in a 3-year orbit around DE0630-18* 

J. Sahlmann^’**, P. F. Lazorenko^, D. Segransan^, E. L. Martin^, M. Mayor^, D. Queloz^’^, and S. Udry^ 

' European Space Agency, European Space Astronomy Centre, P.O. Box 78, Villanueva de la Canada, 28691 Madrid, Spain 
e-mail: johannes. sahlmann@esa. int 

^ Main Astronomical Observatory, National Academy of Sciences of the Ukraine, Zabolotnogo 27, 03680 Kyiv, Ukraine 
^ Observatoire de Geneve, Universite de Geneve, 51 Chemin Des Maillettes, 1290 Versoix, Switzerland 
INTA-CSIC Centro de Astrobiologi'a, 28850 Torrejon de Ardoz, Madrid, Spain 
^ University of Cambridge, Cavendish Laboratory, J J Thomson Avenue, Cambridge, CB3 OHE, UK 

Received 27 January 2015 / Accepted 26 Eebruary 2015 


ABSTRACT 

Using astrometric measurements obtained with the EORS2/VLT camera, we are searching for low-mass companions around 20 nearby 
ultracool dwarfs. With a single-measurement precision of ~0.1 milli-arcseconds, our survey is sensitive to a wide range of companion 
masses from planetary companions to binary systems. Here, we report the discovery and orbit characterisation of a new ultracool 
binary at a distance of 19.5 pc from Earth that is composed of the M8.5-dwarf primary DE0630-18 and a substellar companion. The 
nearly edge-on orbit is moderately eccentric (e = 0.23) with an orbital period of 1120 d, which corresponds to a relative separation in 
semimajor axis of approximately 1.1 AU. We obtained a high-resolution optical spectrum with UVESA^LT and measured the system’s 
heliocentric radial velocity. The spectrum does not exhibit lithium absorption at 670.8 nm, indicating that the system is not extremely 
young. A preliminary estimate of the binary’s physical parameters tells us that it is composed of a primary at the stellar-substellar 
limit and a massive brown-dwarf companion. DE0630-18 is a new very low-mass binary system with a well-characterised orbit. 

Key words. Stars: low-mass - Brown dwarfs - Planetary systems - Binaries: close - Astrometry 


Using astrometry with the FORS2 optical camera (Appen- 
zeller et al. 1998), installed at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO), we are searching for 
planetary companions of 20 southern ultracool dwarfs with spec¬ 
tral types M8-L2. The project is described in Sahlmann et al. 
(2014b), and its first result, the discovery of a low-mass com¬ 
panion to an L dwarf, is reported in Sahlmann et al. (2013) and 
updated by Sahlmann et al. (2015). Details on the astrometric 
reduction methods are given in Lazorenko et al. (2014). 

Here, we report the discovery of the binary nature of DENIS 
J063001.4-184014, hereafter DE0630-18, revealed by ground- 
based astrometric monitoring over more than three years. The 
primary is a very-low-mass star or brown dwarf with an optical 
spectral type of M8.5 (Phan-Bao et al. 2008), and the companion 
responsible for the astrometric orbit is a brown dwarf. 


1. Observations and data reduction 

DE0630-18 was observed as part of our astrometric planet 
search and the FORS2 data were reduced as described in La¬ 
zorenko et al. (2014). The position of the target was repeatedly 
measured relative to reference stars in the 4'x 4' field of view. 
DE0630-18 has an /-band magnitude of 15.7 and we observed 
it on 23 epochs over a timespan of 1209 d between 2010 Decem¬ 
ber 7 and 2014 March 30. On average, we obtained 44 frames per 

* Based on observations made with ESO telescopes at the La Silla 
Paranal Observatory under programme IDs 086.C-0680, 088.C-0679, 
090.C-0786, and092.C-0202. 

** ESA Research Eellow 



Fig. 1. Photocentric orbit of DE0630-18 caused by the gravitational 
pull of the orbiting brown dwarf. Observations and the best-fit model 
are shown as black circles and a grey curve, respectively. Uncertainties 
are smaller than the symbol size. The barycentre and periastron position 
are marked with a cross and an open square, respectively. The orbital 
motion is clockwise; north is up and east is left. 


epoch. In June 2012, it became clear that the standard astromet¬ 
ric model is not sufficient to explain the motion of DE0630-18, 
and we initiated follow-up observations. We followed the same 
strategy and procedures for the adjustment of an additional Kep- 
lerian motion that is described in Sahlmann et al. (2013). A ge¬ 
netic algorithm was used to efficiently explore the large parame¬ 
ter space and to identify the most promising model parameters. 
These provided the starting values for a Markov-Chain Monte 
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Table 1. Orbital parameters of the DE0630-18 system. 
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Notes. Parameter values correspond to the median of the marginal pa¬ 
rameter distributions and uncertainties represent Icr ranges. The proper 
motions are not absolute and were measured relative to the local refer¬ 
ence frame. 

Carlo (MCMC) analysis that yielded the final parameters and 
their confidence intervals and correlation estimates. 

We also observed DE0630-18 on 2013 October 3 (MJD' 
56568.350793) with the red arm of UVES at the VET (Dekker 
et al. 2000) using a 1'.'2 slit width, which provided a resolving 
power of R~33 000, and the Dichroic 2 standard setup centered 
at 760 nm to cover the wavelength range of 565 - 931 nm. The 
exposure time was 2100 s and the observation took place with 
0777 optical seeing at an airmass of 1.10. The spectrum was re¬ 
duced using the ESO pipeline in standard setup. We measured 
the radial velocity on the UVES spectrum using strong atomic 
lines (Rbi and Csi) with the method described in Sahlmann 
et al. (2015) and determined a heliocentric radial velocity of 
-13.0+ 1.1 kms-'. 

2. The photocentre orbit of DE0630-18 

The photocentre orbit of DE0630-18 is shown in Eigs. 1 and 4. 
Compared to the measurement precision, the photocentre semi¬ 
major axis of 23.4 milli-arcseconds (mas) is very large. The orbit 
is seen almost edge-on. Eigure 2 shows the proper and parallac¬ 
tic motion of the DE0630-18 system relative to the background 
stars, which is one order of magnitude larger than the orbital mo¬ 
tion. 

The orbit fit parameters and their confidence intervals are 
reported in Table 1, where Aa^ and Ado relative offsets to 

the target’s position at the reference date taken as the arith¬ 

metic mean of the observation dates, m is the relative parallax, 
Pa* and ps are the proper motions, e is the eccentricity, cj is 

' Modified Julian date (MJD) is barycentric Julian date - 2400000.5. 



Offset in Right Ascension (mas) 

Fig. 2. Proper and parallactic motion of DE0630-18 relative to the 
field of reference stars. The astrometric observations and the model are 
shown as black circles and grey curve, respectively. The black arrow in¬ 
dicates the direction and amplitude of the proper motion over one year. 
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Fig. 3. Marginal distribution of the orbital period (left) and inclination 
(right) for DE0630-18 obtained from 1.5 x 10® MCMC iterations. The 
dashed line indicates the median value and the Y-coordinate units are 
arbitrary and indicate relative occurrence. 

the argument of periastron, P is the orbital period, dRef is the 
mean longitude at time TRef, Q. is the ascending node, i is the or¬ 
bit’s inclination, and a is the semi-major axis of the photocentric 
orbit. The parameters p and d model the differential chromatic 
refraction and Sa and are nuisance parameters. The astromet¬ 
ric data constrains all orbital parameters well, which leads to 
normally-distributed parameters and small fractional uncertain¬ 
ties, see Fig. 3. The parallax correction Am was determined in 
Sahlmann et al. (2014b) and yields the absolute parallax UTabs- 
The time of periastron passage Tq can be retrieved via the mean 
anomaly 

M = T-w = 27r(^-<^o) ^ ro = rRef-P^, (1) 

where t is time and 0o = Tq/P is the phase at periastron. 

On the basis of the FORS2 astrometry, we have accurately 
determined the photocentre orbit of DE0630-18, in particular 
its semimajor axis and period. Additionally, we measured the 
absolute parallax, thus know the system’s distance from Earth. 
Despite the nearly edge-on configuration, the binary is unlikely 
to eclipse because of the large orbital separation between two 
approximately Jupiter-sized bodies. 
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Fig. 4. Photocentre orbit motion of 
DE0630-18 as a function of time. The or¬ 
bital signature in Right Ascension (panel a) 
and declination (panel c) is shown, where 
black symbols show the epoch average val¬ 
ues. Panels b and d show the observed mi¬ 
nus calculated (0-C) residuals of epoch av¬ 
erages. 


3. Preliminary constraints on the individual 
components 

In Sahlmann et al. (2014b), we estimated a mass of Mi - 
0.086 + 0.009 Mq for the primary DE0630-18A. Should the sys¬ 
tem’s photocentre and barycentre coincide, the measured orbit 
corresponds to a mass of M 2 - 0.060 Mq for the companion 
DE0630-18B. However, at this mass ratio q - M 2 /M 1 ^ 0.70, 
the optical light contribution of the secondary can be significant 
and the photocentric and barycentric orbits may be different. 

The difference between a binary’s photocentre and barycen¬ 
tre orbit size is determined by the magnitude difference Am be¬ 
tween the components and their individual masses. The frac¬ 
tional mass / = M 2 /(Mi -h M 2 ) and the fractional luminosity 
P — L 2 l{Li + L 2 ) = (1 -H 10°'*^'”)“*, define the relationship 
between the semimajor axis a of the photocentre orbit and the 
semimajor axis Orei of the relative orbit, where both are measured 
in mas: 

O'= arei(/-jS). (2) 

An independent constraint on the relative semimajor axis is 
given by Kepler’s law 

G(Mi+M2)-4;r2^, (3) 

where G is the gravitational constant, drei is measured in metres 
and P is in seconds. The relation between drei and Orei is given 
by the parallax nJabs- 

Because we have knowledge of a, P, and only, the 
problem is underconstrained. To obtain a preliminary estimate 
of the possible values for individual masses, we therefore have 
to make use of theoretical mass-luminosity-age relationships. 
Eor a range of companion masses M 2 and a constant primary 
mass Ml, we used the BT-Settl (Allard et al. 2012) and DUSTY 
(Chabrier et al. 2000) models to obtain the corresponding mag¬ 
nitude difference Am/ in the /-band for ages >1 Gyr. Because we 
have measured a, Eq. (2) yields the estimate for the relative 
semimajor axis. On the other hand, we know the orbital period 
and the parallax, thus Eq. (3) gives us a second estimate a"j of 
the relative semimajor axis. The only possible values of M 2 are 
the ones where the identity a"j = a'^j is fulfilled. 

In Eig. 5, we show the relations between a"j and a'^j for dif¬ 
ferent system ages. The dashed grey line indicates equality and 
isolates two values of M 2 that are allowed for every age. The 
main drawback of this method is that it relies on models that 
may not be well calibrated in the respective mass and age range. 
However, it allows us to draw first preliminary conclusions on 
the system parameters. 



Fig. 5. Relative semimajor axis of DE0630-18 obtained from two es¬ 
timators using the mass-luminosity-age relationships according to BT- 
Settl (black curves) for ages of 1 Gyr (dashed), 3 Gyr (solid), and 7 Gyr 
(dotted) and DUSTY (thick grey curves) for ages of 1 Gyr (dashed) and 
5 Gyr (solid). The secondary mass increases from 0.050 Mq to 0.076 Mq 
upwards along every curve. Allowed values are marked with open and 
filled circles. 


According to Fig. 5, the system has to be older than 1 Gyr be¬ 
cause the requirement a"j = cannot be met at this age. This 
is compatible with the non-detection of Lit absorption in our 
UVES spectrum of DE0630-18. Because objects with masses 
>0.06 Mq deplete their lithium content within the first hundred 
million years of existence, the lithium test (e.g. Magazzu et al. 
1993) yields a constraint on mass and/or age of ultracool dwarfs. 
For DE0630-18, the Li test is negative, thus DE0630-18 is not 
very young, in agreement with the constraint from the orbit mod¬ 
elling. 

Possible secondary masses are found at ages of 3, 5, and 7 
Gyr, where the secondary can have a mass of M 2 - 0.075 Mq 
or of M 2 - 0.060 Mq. These two solutions correspond to sce¬ 
narios where the companion is very massive and luminous, thus 
contributes significantly to the photocentre-to-barycentre shift, 
and where the secondary is less massive, leading to a smaller 
orbit and smaller photocentre-to-barycentre shift, respectively. 
Table 2 lists the model-dependent values of secondary mass, rel¬ 
ative separation, magnitude difference in red and infrared bands, 
and the difference between the photocentre orbit size and the pri¬ 
mary’s barycentric orbit size (a-ai) for all solutions. The binary 
separation is in the ~57-59 mas range and the magnitude differ¬ 
ences range between 2.9 and 13.3 in the /-band and between 1.1 
and 5.3 in the infrared K-band. 
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Table 2. Theoretical constraints on secondary masses and magnitudes 
according to BT-Settl and DUSTY models. 
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59 
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5 
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57 
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4. Discussion 

Astrometric monitoring of the ultracool dwarf DE0630-18 led 
to the discovery of its binary nature and allowed us to determine 
all orbital parameters with high precision. Because the photo¬ 
centric motion in /-band may be diluted by the companion’s 
emission, additional constraints are required to characterise the 
individual binary components. Theoretical models of substellar 
evolution point towards system age older than 1 Gyr and a pri¬ 
mary mass of Ml = 0.086 + 0.009 Mq. The secondary mass 
has two allowed modes at ~0.06 Mq and ~0.075 Mq, which in 
both cases indicates a substellar nature, i.e. the companion is 
a brown dwarf. From our astrometric and spectroscopic obser¬ 
vations alone, we cannot distinguish between these two modes. 
However, we note that DE0630-18 is not overluminous in the 
Wide-held Infrared Survey Explorer (WISE; Wright et al. 2010) 
bands Wl, W2 (shown in Fig. 6), and W3 compared to the popu¬ 
lation of ultracool dwarfs, which points towards the lower mass 
for the companion. In Fig. 6, the juvenile binary DE0823-49 
(LI.5, Sahlmann et al. 2015) appears slightly overluminous. The 
seemingly overluminous object is DE1733-16 (El), which is 
probably caused by blending with a background source in this 
extremely dense starheld. 



Fig. 6. Absolute magnitude in the W2-band (~4.6 pm) as a function of 
spectral type for M6-L5 dwarfs in the database of ultracool parallaxes 
(Dupuy & Liu 2012) (grey symbols) and for our survey sample with 
parallaxes from Sahlmann et al. (2014b) (black symbols). Magnitude 
uncertainties of the latter are smaller than the symbol size. Tight bina¬ 
ries in our sample are shown with open circles. DE0630-18 (M8.5) is 
seen to have a magnitude similar to other ultracool dwarfs. 


Figure 7 shows the expected radial velocity curve in this 
scenario. Our UVES measurement was taken close to a time 
of zero orbital radial velocity, which allows us to determine 
the heliocentric systemic velocity as either -12.2 + 1.1 kms^' 
or -13.8 + 1.1 kms“* (for the opposite orbit orientation with 
i' = i -I- 180°). Clearly, the radial velocity variations along the 
orbit are detectable with follow-up measurements using UVES 
or other spectrographs. 



Year 


Fig. 7. Estimated radial velocity curve of DE0630-18A (Mi = 
0.086 Mg, M2 = 0.061 Mg) as a function of time. The UVES measure¬ 
ment is shows with a solid circle and the systemic velocity is indicated 
by the horizontal dashed line. The dash-dotted line corresponds to the 
alternative radial velocity curve with i’ = i + 180°. 

A more detailed characterisation of the binary components of 
DE0630-18 requires additional observations. The relative prox¬ 
imity (19.5 pc) and brightness (;«/ ^ 15.7, nij ^ 11.3) of this 
system will facilitate these and possibly lead to new insights into 
the properties of ultracool dwarfs. 

The case of DE0630-18 illustrates the potential difficulties 
when trying to characterise ultracool binary stars from the pho¬ 
tocentric orbit alone, which will also arise for the hundreds of 
orbits (Sahlmann et al. 2014a) expected from the astrometric sur¬ 
vey of the Gaia mission. 
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